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This is an extended version of the note taken by the first author (W.-G.P.) on a lecture 
given by the second author (M.R.) as a first part of the series on "Hadronic Matter Under 
Extreme Conditions," the principal theme of the WCU-Hanyang Program. It covers the 
attempts to go in a framework anchored on effective field theory of QCD from zero 
density to the nuclear matter density and slightly beyond, with the ultimate goal of 
arriving at the density relevant to compact stars, including chiral phase transition and 
quark matter. The focus is on the conceptual aspects rather than detailed "fitting" of 
the data on the kinds of physics that are being addressed to in radioactive-ion-beam 
machines in operation as well as in project (such as 'KoRIA' in Korea) and will be 
explored at such forthcoming accelerators as FAIR/GSI. The approach presented here is 
basically different from the standard ones found in the literature in that the notion of 
hidden local symmetry - which underlies the chiral symmetry of the strong interactions 
— and its generalization to dual gravity description involving infinite tower of hidden 
gauge fields are closely relied on. 
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1. Introduction 

What we would like to do in this lecture is to figure out what happens in dense 
matter at zero temperature expected in compact stars. This lecture is about what 
we understand up to the present time, and what we would like to understand in the 
future. 

We start with a cartoon panorama. 

*The World Class University (WCU) Lecture given in the Fall 2009. 

1 



16, 2010 7:51 WSPC/INSTRUCTION FILE 



WGP-MR-arXiv 



2 Won-Gi Paeng and Mannque Rho 



T(emperature) 



Quark/gluon plasma (?) 



Early Universe 




RHIC etc 



> 

Compact stars & black holes 



n (density) 



Fig. 1. The phase diagram of hadronic matter in temperature and density that is being studied. 
The solid line describes experimental effort, e.g., SPS/CERN, RHIC and the dotted line lattice 
QCD. 

We are often shown, very roughly, a phase diagram like Figure 1, which is given 
in terms of density vs. temperature. When you go along the temperature axis (and 
at zero density), you are probing the Early Universe. Along the temperature axis at 
near zero density, nothing drastic happens. There is no real phase transition of the 
kind naively predicted and widely heralded by theorists in the past. This feature 
is confirmed by lattice QCD calculation. Experiments are probing a wide range 
of temperature and density (indicated by solid line in Figure 1). We don't know 
actually where what one would call genuine "quark gluon plasma" (QGP) is. But 
one expects that QGP should be there somewhere in the figure. In a phase where 
there are weakly or non-interacting quarks and gluons with asymptotic freedom, you 
can calculate things in perturbation theory. So thinking that one will get quickly to 
QGP, people have done such calculations, a huge number of them. 

Many people were dead wrong. The naive QGP is not visible yet. It's a lot more 
subtle. Whatever it is, it is strongly coupled up to quite high temperature and does 
not even require QCD proper to describe it. Some sort of universality involving 
conformality is in action, blissfully ignoring whether certain symmetries such as 
supersymmetry are there or not. This seems to mean that one is not really learning 
anything significant here. 

Furthermore, whatever physics it is that you are describing in terms of quarks 
and gluons of QCD, at low energy regime, you can also describe in terms of hadrons 
and other degrees of freedoms which are actually measured in experiments. Lots of 
people are calculating what happens in low-density non-zero temperature regime 
which are being and will be explored by experiments, at e.g., SPS/CERN, RHIC, 
ALICE/LHC etc. Here lattice plays a crucial role in guiding theorists. 

But lattice has tremendous difficulty handling density and in particular, the 
high density appropriate for compact stars seems out of reach of lattice techniques. 
There are no experiments to guide theorists beyond the normal nuclear matter 
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density uq - and the experimental facilities for this regime are yet to come, so there 
are no trustful model approaches either. We did not put in the figure anything novel 
dealing with density. You see why. It is a fair picture of the current status. There 
is nothing much known there. Dense matter is a totally virgin field to explore and 
as indicated on Fig. 2, this is the WCU objective. The only thing we are sure of is 
the nuclear matter at n = uq = 0.16 fm -3 . There is also a large effort at various 
laboratories in the world, e.g., GANIL and RIB machines (e.g., also forthcoming 
KoRIA?), to unraveling the gas- liquid phase transition at lower density. 

2. Phenomenology with Increasing Density 



"Strange" things start 
happening here 

4- 

Fermi liquid 

-[ 1 — n - J 1 n = Number N/Volume 



gas-liquid n ° 

WCU 



t- 



Fig. 2. First there is the gas-liquid phase transition point, then the (Fermi) liquid phase. Nuclear 
physics up to the nuclear matter density no = 0.16fm is very well established. From no up to 
the chiral restoration density n x is the density domain for the WCU-Hanyang Program. 



2.1. Nuclear physics up to ~ 3no 

Without much ado, let's jump into a specific case raised by Bao-An Li's talks given 
here. This case illustrates all the interesting features we are particularly interested 
in. Figure 3A shows the behavior of nuclear symmetry energy - which figures im- 
portantly in neutron-rich nuclei, compact stars etc. - at density above n in various 
different models QEEl) . It shows how widely diverging are the predictions for the 
symmetry energy given by various models that all fit nuclear matter properties up 
to, and perhaps slightly above, nuclear matter density. The divergence simply shows 
how little constrainable the symmetry energy is both by theory and by experiment 
at higher density. The situation is not so bad for the energy for symmetric (N = Z) 
part (although beyond uq there is certain spread in prediction), so we will focus 
on the symmetry energy that enters where there is neutron excess which is the 
tendency in Nature, both in laboratories and in compact stars. 

In most of models in the literature, the symmetry energy increases steadily with 
increasing density. But the one that fits the available ir~/ir + ratio measured by 
FOPI/GSI and analyzed by Bao-An Li El - there is only one data at the moment 
- differs widely from the rest (call it "canonical" ) in that it increases up to certain 
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Fig. 3. Diverse symmetry energies predicted by different models (A) which are caricatured by 
the simplified one (B).HS1. 



density, say, around ~ 1.3rio and then turns over and decreases steadily, going to 
zero at n ~ 3nrj. Let us call it "supersoft" following Bao-An Li. For convenience, we 
just replace this complicated picture by the simplified one Fig. 3B. We will address 
the problem referring to Figure 3B. 

Suppose Li's symmetry energy S a —i (called "supersoft") is what's given by 
Nature. There are then several serious consequences: 

(1) With the standard description of kaon- nuclear interactions ^ - largely based 
on chiral perturbation theory - which predicts that an anti-kaon mass drops 
continuously in compact star matter as density increases, there cannot be kaon 
condensation at any reasonably low enough density to be relevant to neutron- 
star physics. 

(2) Unless one invokes non-Newtonian gravity, there would be no stable neutron 
stars of the range of masses and radii observed so far, clearly at odds with 
nature^. 

(3) Then all the "exciting things" like the Bethe-Brown maximum mass of neutron 
stars, the maximization of black holes in the Universe and the cosmological 
natural selection discussed in ^ will be nothing but red herrings. 

These issues will pose great challenges to nuclear physicists. The "standard" chiral 
perturbation approach typically gives the kaon mass behaving like what's shown as 
"%PT" in Figure 4. The kaon with that sort of dropping mass cannot possibly con- 
dense in compact-star matter, since the electron chemical potential [i e is not high 
enough to make the "electron decay" e — > K~ + v. Such a "supersoft" symmetry 
energy will typically make a havoc to neutron stars as one can verify by sticking it 
into the TOV equation. As Li suggests, one can make drastic measures to "save" the 
soft symmetry energy. One way-out suggested by Li is to introduce non-Newtonian 
gravity without disturbing what's established with standard gravity theory. This 
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spells like a disaster but it could also be an exciting opportunity for the futureQ 
Because when something goes wrong, there is a ground for a breakthrough: One 
now has to find out whether the S a =i is what Nature actually says. Experimen- 
tally, a part of the symmetry energy - the lower density part - could be measured at 
radioactive ion beam (RIB) machines in construction all over the world, including 
the one Koreans are going to build, called "KoRIA," and with precision measure- 
ments, one could check whether the FOPI experiment is actually constraining the 
symmetry energy that drastically. There are also big theoretical issues, the solu- 
tions of which could provide ways out of the disaster. We will discuss below two 
related possibilities based on BR scaling^ anc j three-body (and many-body) forces 
that would solve the problem even if the supersoft symmetry were correct up to the 
measured density. 




Fig. 4. The electron chemical potential /i e which is directly related to the nuclear symmetry 
energy corresponding to "supersoft," the "standard" kaon mass labeled %PT predicted by the 
class of models belonging to chiral perturbation theory and the kaon mass predicted in dense 
half-skyrmion matter. The "standard" kaon mass does not meet with the chemical potential given 
by "supersoft," so there will no kaon condensation. However the half-skyrmion phase makes the 
kaon mass drop sufficiently fast that it could trigger kaon condensation. 



Suppose kaons do condense in compact stars at a density near 3no- Below we will 
discuss how this can happen because of some highly non-perturbative mechanism 
recently discovered. If this happens, it has an interesting implication on laboratory 
experiments. 

Kaon condensation can take place because there is effectively strong kaon-nuclear 
attraction. Toshi Yamazaki, who is a well-known Japanese nuclear experimentalist 
whom the lecturer (MR) admires (one of those rare nuclear experimentalists who 
do not take seriously what theorists say), predicted that if the anti-kaon nuclear 



a Indeed, it is a fashionable idea that Newtonian gravity arises as an "emergent" phenomenon 
through a holographic scenario. This possibility offers an avenue for a conceptual breakthrough 
for string theorjEl as we n as f or loop quantum gravity theorjEl. 
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potential were as deep as 

Vk-jv(no) £-200MeV (1) 

then there could be highly bound kaon-nuclear systems with an average density 
much higher than that of nuclear matter. His argument^] as j s being debated, is 
not a "kosher," as some of the chiral perturbation theorists who do extensive chiral 
perturbation calculations with various approximations do not agree with his predic- 
tion. But one should recognize that chiral perturbation theory, which is a QCD for 
certain processes (involving pseudo-Goldstone boson) in some low-energy kinematic 
domain^, may not be able to access high density regime. In fact, one doubts that 
what comes out of chiral perturbation theory at low order and involving a limited 
number of degrees of freedom could be sufficiently reliable. As we will discuss below, 
there can be certain nonperturbative mechanisms that cannot be accessed by a stan- 
dard chiral perturbation technique, namely, a phase change of the baryonic matter 
into a half-skyrmion or half-instanton matter with chiral symmetry "restored" in 
the sense that the quark condensate vanishes but with pions fluctuating on top of 
the modified ground state (or "vacuum") and with quarks/gluons still confined in 
hadrons, contrary to the standard scenario. In such a phase, kaons can become very 
light and could condense at low density. 

2.2. Half- sky rmions and nuclear physics at ~ 3no 

As was explained in previous lectures, in the absence of first-principle calculations 
anchored on QCD in approaching dense matter, the best one can do is to start with 

II Q| 

hidden local symmetry Lagrangian^ 21 . It would be much more realistic to do it with 
the 5D Yang-Mills action supplemented with the Chern-Simons term, either decon- 
structed from low-energy or reduced from string theory from high energy. We have 
had a number of lectures on this latter matter in this WCU series. Unfortunately 
we cannot yet handle the 5D theory for technical reasons we have heard quite a 
bit about, such as 1/N C , 1/A corrections etc., so let us do the 4D theory, i.e., the 
Harada-Yamawaki (HY) HLS theory. 

Given the HY HLS Lagrangian, we need to get baryons, and this we can do 
via topology. The baryons are of course needed for sourcing density. Alternatively 
one may resort to Weinberg's "folk theorem" and start with explicit baryon fields 
and HLS meson fields put in a way fully consistent with chiral symmetry to build 
dense matter. This approach will be taken up later. What we will do here is to work 
with these skyrmions generated from the HLS theory. At the moment, there is no 
continuum treatment using the skyrmions combined with HLS mesons 0. One can, 
however, put skyrmions on crystal lattice and squeeze the system to generate high 

b There are attempts to approach dense matter analytically in holographic QCD. In terms of multi- 
instantons as multi-baryons, one can address such questions as many-body forces at short distance 
relevant to high density, e.g., Hashimoto et alEl. This approach is being developed and is not yet 
in a workable form. 
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Fig. 5. Phase diagram near chiral restoration at high density. The relevant densities are in ar- 
bitrary units as they depend on the parameters of the Lagrangian which are to be fixed by the 
physics involved. See Figure 6. 

density. This matter was discussed by Byung-Yoon Park in his colloquium here and 
also extensively reviewed else where, s o we won't repeat it. What comes out of such 
treatments is shown in Figure Jp3HSI At some density (see later for numbers), there 
is a phase change from a skyrmion matter - which is a nucleonic matter - to a half- 
skyrmion matter with vanishing quark condensate and non-vanishing pion decay 
constant. This is a phase in which chiral symmetry is apparently restored (with 
the vanishing order parameter (qq) = 0), but there are hadrons, with presumably 
parity-doubled and massive baryons and mesons, and the Goldstone pions. El Note 
that this is a novel phase nobody else has been discussing up to now although some 
people are arriving at something similar using large N c arguments, which we think 
is more of a conjecture than anything else at the moment. 

Now you have a dense matter with a particular phase structure which is novel. 
What happens when you inject a K~ into the system? This sort of experiments, if 
not done yet, will be done at the J- PARC and surely at the FAIR/GSI. What hap- 
pens to an ani-kaon when embedded in the half-skyrmion matter has been worked 
out. Although the final result appeared after the lecture, let us quote the result here 
since it is rather striking^!. It is given in Figure 6 for the dilaton mass 720 MeV. 
We do not really know what the dilaton mass we need for our study is. The story of 
scalar mesons in low-energy strong interaction physics is still very murky and there 
is no simple understanding. There can be a highly complex interplay between quark 

c It's possible that the 't Hooft anomaly matching is satisfied by the masslcss pion. 
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degrees of freedom and gluonic degrees of freedom. What is given in Figure 6 is for 
the dilaton mass that is used in the rather successful mean-field calculation of nu- 
clear matter performed by Chaejun Songl3 using an HLS Lagrangian in which the 
vector manifestation, i.e, BR scaling, is implemented in consistency with Landau 
Fermi liquid th eory . 

It turns outP^ that the threshold density for the half-skyrmion phase rii/2 ~ 
1.3nn is more or less independent of the dilaton mass. It is controlled by the pa- 
rameters of the skyrmion Lagrangian itself. What the dialton mass does is to shift 
the chiral transition density n x : the higher the mass is, the higher n x mwill be. 
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f*,*0 



Kaons will condense 
in compact stars? 



, V2 skyrmion 



Nuclear matter 
J-Parc 



n/n_ 



Kaons condense 
in nuclear matter 



Fig. 6. Kaon spectrum in the half-skrymion matter for the dilaton mass of 720 MeV. 



The important feature to notice is that in a half-skyrmion phase, the kaon mass 
drops considerably more rapidly than in the usual skyrmion matter, going to zero 
when the condensate of the soft dilaton - which coincides with the chiral conden- 
sate within the approximations taken - vanishes. There seems to be, in the half- 
skyrmion matter, the "mysterious" agent in Yamazaki's picture that shrinks the 
matter, a mechanism that is not readily available in ordinary nuclear interactions. 
Interestingly, kaon condensation can take place even in a symmetric nuclear matter 
when the scale symmetry breaking tied to chiral symmetry breaking is restored a 
la Freund and NambvP^, implying that kaon condensation is equivalent to strange 
quark matter. If one takes the canonical symmetry energy of El, then kaon conden- 
sation will take place at n ~ 3nrj and things will happen as expected. In fact, the 
results obtained in this calculation are s urprisingly close to what was obtained in 
the HLS/VM scenario described in EE. What will happen with Li's "supersoft" 
symmetry energy is not clear. It's a nice problem to work out. 
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3. Nuclear Matter in "Double Decimation" 

Going to higher densities in nuclear systems involves a multitude of scales. When 
there are a multitude of scales involved in many-body processes as in the case 
of dense-matter physics, one needs to figure out which is the best starting point 
for a given process. It is clearly senseless to start with physics at nuclear matter 
density and hope to get, within the same framework, to super-high density where 
perturbative QCD is applicable, such as for instance the color-flavor locking phase. 
There is a good example at high temperature as studied at RHIC. The finding - 
which is a total surprise to many people working in heavy-ion physics - was that a 
conformally invariant supersymmetric field theory which has nothing to do with low- 
temperature world can describe well what's happening at high temperature where 
they all thought weakly interacting quark-gluon plasma would show up. It appears 
that whatever symmetries and dynamics that are essential at high temperature for 
the processes involved there have little to do with what happens at low temperature 
which can be looked at by lattice QCD. It may be that similar things happen in 
density. 

One thing is certain, and that is that we should be able to go with confidence 
from zero density up to near the nuclear matter density. Let's see what we can say. 
Most of the workers in the field start with a Lagrangian determined in matter-free 
space at zero temperature. Can such a theory go beyond the nuclear matter density 
and go up to the chiral transition point? 

We know that the ordinary nuclear matter we are familiar with is a fixed point 
of an effective field theory (EFT), namely, Landau Fermi liquid fixed point. Ignor- 
ing the possibilities of meson condensation, half-skyrmion (or half-instanton) phase 
etc, there is next the chiral transition point which in HLS theory is the "vector 
manifestation (VM) fixed point." The question is: Can one go from one fixed point 
to another in a systematic way and if so, how? 

It is clear that by just sticking around the matter-free space and insisting on 
calculating more and more terms - which is what many nuclear theorists are doing 
- one will not get anywhere. One has to find how to proceed in a consistent and 
efficient way. We don't know how to do this all the way to the critical density for 
chiral restoration - and that's why we have been playing around with skyrmion 
crystals, vector manifestations etc, but we can do it up to, and slightly above, the 
nuclear matter density. And this is what we will do in this part. How to go from 
nuclear matter above that density is something that our WCU program will have 
to figure out. 

We will get to nuclear matter by what Gerry Brown and MR called "double 
decimation." For this the first thing is Weinberg's "folk theorem" which was restated 
this year in 

3.1. Weinberg "folk theorem" 

The theorem goes as follows for the case we are concerned with: 
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"If we write down the most general possible Lagrangian, including all terms 
consistent with chiral symmetry principle, and calculates matrix elements with the 
Lagrangian in any given order of perturbation theory, the result will simply be the 
most general S-matrix consistent with perturbative unitarity, analyticity, cluster 
decomposition and the assumed symmetry properties" 

We will now sketch how this theorem guides us to understand physics of nuclei 
and nuclear mater. 

3.2. Double decimation 

If we are interested in what's happening in the vicinity of the nuclear matter satu- 
ration point, then we can simply start at that point taken as a Fermi liquid fixed 
point as we will do below. One can write down an effective Lagrangian with correct 
symmetries, principally chiral symmetry, and calculate weak fluctuations around the 
Landau Fermi liquid fixed point. This can be phrased using an effective Lagrangian 
patterned after Walecka's mean field theory. One starts with a suitable HLS La- 
grangian with both fermionic and bosonic degrees of freedom that are relevant and 
work in the mean field. This is shown to be equivalent to doing Fermi liquid theory 
as we will explain further below. Of course it's not obvious how to go beyond the 
vicinity of the normal nuclear matter, but the physics just below and just above nu- 
clear matter density should be describable as Walecka's mean-field theory does. One 
may extrapolate beyond the nuclear matter density in some sense if one endows the 
parameters of the theory with density dependence. This renders the theory a sort of 
density functional theory of Kohn-Sham type in field theory language. This is the 
strategy of BR scaling - something which has to do with chiral symmetry - making 
contact with standard nuclear physics which does not have obvious connection to 
chiral symmetry. It turns out that this works very well. Bingo for BR scaling! 

The question that can be answered is this: Can one go from a Lagrangian defined 
in the matter-free vacuum to the Fermi liquid fixed point? The answer is yes: It 
involves the double decimation. 

Suppose we have an effective Lagrangian C e f / which has all relevant degrees of 
freedom and symmetries required by QCD defined with a cut-off at A. In practice, we 
are thinking of chiral Lagrangian for which the cut-off is A = A x w 4^/^ ~ 1 GeV. 
One usually assumes that in dealing with many-baryon systems, there is a Fermi sea 
defined by the Fermi momentum hp. In principle, the Fermi sea should come out of 
the theory, as most likely a quantum critical phenomenon. The only discussion on 
this matter in co nne ction with nuclei and nuclear matter as far as we know was made 

lOI I 

by Lynn in 1993 where he suggested that a drop of nuclear matter, i.e., nucleus, 
could arise as a non-topological soliton, a chiral liquid from the chiral Lagrangian 
with which the effective potential is calculated in chiral perturbation theory. In 
practice, one cannot do this to all orders since an effective theory defined with a 
cut-off would require an infinite number of parameters for an all-order calculation. 
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Even when drastically truncated, the theory could not quite arrive at a form that 
resembles the requisite Fermi liquid, and no further work along the line initiated by 
Lynn has been done since. 



vacuum (Q,m JI )/A x 



Fig. 7. Multiple scales and expansion parameters. Here Q is the typical momentum/energy scale 
probed, m n is the pion mass due to the explicit breaking of chiral symmetry by the quark mass. 
Note that up to A p , the Fermi momentum appears as a "small" expansion parameter but near 
the Fermi surface, it appears as a "large" N in 1/N expansion in the Fermi-liquid fixed point 
Lagrangian. 

A program that has made a big progress is that of Gerry Brown and his collabo- 
rators. They assume that there is a Fermi sea and derive nuclear matter in two steps, 
in what Gerry and MR called "double decimation." This is indicated in Figure 7 
with the pertinent expansion parameters. The first step is from zero energy scale 
(A = 0) roughly to the scale up to which nucleon-nucleon (NN) elastic scattering is 
measured and the second is from the free-space interactions to the Fermi sea scale 
with many-body interactions. Here one is doing an effective field theory, making the 
expansion in terms of some small expansion parameters. Now there is an intricacy 
in the way the expansion is to be done in the two steps. The first resorts typically to 
chiral perturbation theory where the expansion parameters are Q/A and m T /A with 
Q being the small probe momentum, m v small chiral symmetry breaking and A the 
chiral scale ~ 47r/ 7r . In this step, the Fermi momentum kp is considered as O(Q). 
However in the second step where one will ultimately have to perturb around the 
Fermi surface, kp comes as 1/kp, counted as 0(1/ N) where N is to be considered 
as large. This is indicated in Figure 7. 

Let us describe these two steps in some detail. 

3.3. Vlowk 

Nuclear potential is not a measurable quantity; it is however a convenient inter- 
mediate tool to get to measurable quantities. It's a nice way of translating the 
large amount of NN scattering data given in terms of phase shifts, accumulated by 
experimentalists over the years, into physical observables in many-body systems. 
Assuming that two-body interactions are dominant in many-nucleon systems, the 
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potential extracted from scattering data plays a key role in calculating the proper- 
ties of nuclei and nuclear matter. We will come to n-body forces for n > 2 below. 
Here we will continue ignoring n > 2. 

There are in the literature many such "realistic potentials" constructed from 
the given experimental phase-shift data, more or less guided by certain theoretical 
framework, such as meson exchanges, dispersion relations etc. The examples are 
Bonn A, CD-Bonn, Paris, Nijmegen I and II, Argonne vl8, chiral perturbation, 
Idaho A etc. They are all fine-tuned with a certain number of free parameters El to 
fit the elastic scattering data, typically up to the lab momentum p max ~ 300 — 350 
MeV. The fit is achieved with x 2 /dof « 1, some better than others. One thing 
common among all these "realistic potentials" is that the long-range part has to be 
consistent with one-pion exchange. However the short-range parts can be drastically 
different from one potential to another. If one were to naively use these potentials in 
perturbation expansion to compute the properties of nuclei and nuclear matter, the 
short-distance components of the potential would come in very importantly, which 
account for the tremendous number of papers written on now to do many-body 
calculations in nuclear physics. 

An adroit use of Vi ow k resolves this conundrum. How this comes about is the 
story of Wilsonian renormalization group (RG) technique. 

An effective theory is defined with a cutoff A above which the theory is not 
valid. This is clearly defined in chiral perturbation approaches but with other phe- 
nomenological potentials, this is not the case. To be as precise as possible, we will 
put the arguments in the framework of effective field theory. 

Effective field theories dealing with chiral dynamics have the cutoff at the chiral 
scale ~ 47r/ ff . Suppose we have a two-body potential Vnn determined by fitting 
elastic scattering data up to the momentum p m ax- The NN scattering amplitude A 
is given by the infinite sum 

/t\A\A\ if\\r M\ _l V' (f\VNN\n)(n\V NN \i) 

(f\A\i) = {f\V NN \i) + E- -E ( ' 

n— 1 

where the ellipsis stands for the infinite higher order terms. Here the intermediate 
sum is written to go up to infinity but it should actually be limited to the cutoff 
scale, in principle well-defined in EFT, but in phenomenological potentials, the 
cutoff is done with form factors etc. which cannot be given a precise meaning. Now 
the maximum probe momentum is p m ax up to which precise data are available, so 
for observables probed at p <C Pmax, it makes sense to lower the cutoff from A x 
(or oo) to Ann ~ Pmax- One can rewrite ([2]), defining an effective low-momentum 



d The potentials derived from chiral perturbation theory at low orders have fewer parameters, but 
going to higher order would require an unmanageable number of parameters. 
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potential Vi ow k which results from integrating out high energy mode, as0 

/f\AU\ //it/ i-\ _l (f\ v iowk\n)(n\Vi owk \i) , . 

(f\A\t) = {f\Vt owk \t)+2_j e - E ™ 

n— 1 n 

Note that Ann is arbitrary, so physics should not depend on where Ann is put. 
This gives the renormalization group equation (RGE) 

ffl4!>=o. (4) 

How the effective potential runs with the scale Ann is given by 

d 



,. -Viowk = /3([Viowk], Ann)- (5) 

CIAnN 

The P function is then calculated order by order in an EFT scheme, e.g., chiral 
perturbation. Note that when the (3 function vanishes, one has a fixed point. We 
will have this fixed point in the second decimation. 

This is all fine in principle. But in practice, one is limited by certain calculational 
techniques; the /3 function cannot be calculated to all orders, so it's not simple to 
satisfy Eq. (j4j. The strategy is then to search for Ann where the RGE is best 
satisfied. One finds that this comes out well for Ajvtv ~ 2.1 fm _1 ss 3m„-. 

Given a potential Vjvjv, calculating the corresponding Vi ow k cannot always be 
done analytically. In all cases except in special approximations, it can only be done 
numerically. It is instructive to construct a schematic model for the potential of 
different forms that can be handled analytically which can be fit to experimental 
data in a p articular channel. It turns out that this can be done for the So phase 
shift. In 1221, two separable potentials called I and II were constructed to fit via Eq. ([2]) 
the 1 5'o phase shift. While the Vnn' s so obtained give quite different diagonal matrix 
elements as shown in Fig. 8 (left panel), the Vi ow k's resulting from the decimation 
give identical matrix elements as shown in Fig. 8 (right panel). Also shown there 
are the results obtained with the realistic potentials, which all converge on the same 
curve if one takes Ann = 2.1 fm _1 



3.4. Landau Fermi liquid 

Given the Vi ow k determined in matter-free space from NN scattering, how would 
one go about describing many-nucleon systems, i.e., nuclei and nuclear matter? In 
medium, Vi ow k will get modified by change in the vacuum structure, i.e., by BR 
scaling, and/or rcnormalized by many-body processes. The latter could be calcu- 
lated in a systematic many-body theory. But there also can be intrinsic many-body 



c This is the first decimation in the the double decimation we are doing. 

f For different values of Ajvjv, not all potentials agree in Vi ow k m various channels. It's not clear 
why Ajviv = 2.1 fm —1 gives a better convergence than other values that are somewhat higher or 
somewhat lower. It could have something to do with various other degrees of freedom starting to 
come into the picture. 
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The diagonal momentum-space matrix elements for separable potentials I and II. 



forces, not accounted for in many-body scatterings. For instance in chiral perturba- 
tion theory, there are three-body forces that appear at 0((Q/A x ) 3 ), Figure 9. We 
will not have anything to say how three-body and higher-body forces should figure 
in the RG procedure. There is very little done in that direction. Let us proceed with 
the two-body Vi ow k alone. 

Consider describing many-body systems in which nucleons are interacting via 
Vi ow k in a Fermi sea. We don't know how the Fermi surface arises. We will simply 
assume that there is a Fermi sea formed somehow in the theory. We are interested 
in computing the Bethe-Salpeter equation for fully irreducible vertex V in terms 
of particle-hole irreducible vertex V given in Figure 10. In calculating the integral 
equation, the momentum cutoff is set to Atvat = 2.1 fm _1 . The Landau Fermi liquid 
theory is given by the vertices in the so-called "fc limit" and "w limit" : - — - — > and 
where K = (lj,K)^, 



I A' | 



r w (fci,fc2) = Hm lim T(k 1 ,k 2 ;K) 



w—>0 



A'|->0 



T k (ki,k 2 )= lim limT(k 1 ,k 2 ;K) 



The forward scattering amplitude T k is given by 



Fa8,-ys(kli ^2) — 



^N Z 2 J dO.r^fei, q)T k ^ s {q, k 2 ). 



What is known as "Landau / function" is 

f(k 1 ,k 2 ) = z 2 r u (k 1 ,k 2 ) 

and "Landau a function" - which is forward scattering amplitude 

a(k u k 2 ) = Z 2 T k {k 1 ,k 2 ). 



(6) 

(7) 
(8) 
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Ci,C 3 ,C 4 C'D C E 

Fig. 9. Three-body forces that contribute in chiral perturbation theory at 0((Q/A x ) 3 ). The 
constants Ci,3.4 are determined from 7rN scattering and Cjj can be determined via PCAC from 
the axial current matrix element of two-body systems. The constant Ce is the only undetermined 
parameter in free space. 



In nuclear physics, both spin and isospin figure, so the Landau / function has the 
components (in condensed matter, the isospin terms are absent) 

f(h,k 2 ) = ^ r [F(k 1 ,k2)+F'(k 1 ,k2)T 1 -T 2 + G(k 1 ,k 2 )(7l-(J2+G'(k U k 2 )T 1 -T 2 Cr 1 -Cr2}. 

There are also tensor terms but we will not write them down. 

Given these Landau functions, a small change in the single particle distribution 
function 5n(k) induces the change in energy of the form 

SE = ^ ] 5n{k) + f{ki,k 2 )5n(ki)5n(k 2 ) + 0((5n(k)) 3 ) (9) 

k fci,/C2 

From this one can calculate all the bulk properties of nuclear matter. For instance, 
some relevant physical observables are 

m* F x , . 

— = 1 + ^ io 

TO 6 
h 2 k 2 

+ (") 

K= M_kj L{1 + Fo) (12) 

TO* 

6gi = IH^Il (is) 
y 61 + F1/3 v ' 

where F(k, k') = F1P1 (cos 9). They correspond, respectively, to the effective mass 
of the quasiparticle, symmetry energy, compression modulus and anomalous orbital 
gyromagnetic ratio. 

These Landau / functions have been computed by s umm ing Feynman diagrams 
using the technique of induced interaction with Vi ow k in^l What one finds is quite 
interesting and important. 

When calculated with the Viowk determined in matter-free space (elastic scat- 
tering) and with no 3-body forces, the results (given in the first column of (fl4|)) are 
not quite satisfactory. One notes that the discrepancy in Sgi is particularly serious. 
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Fig. 10. Bethe-Salpeter equation 

3.5. Enter BR scaling and/or three-body forces 

There are two possible ways to improve these results l^sil One is to incorporate BR 
scaling in Vi ow k which is to endow it with density dependence, i.e., Vi ow k{n). The 
way to do this is to consider one-boson exchange potentials - such as Nijmegen 
potential I and II and CD-Bonn potential - and endow BR scaling to the mass of 
the exchanged particles. In doing this, one assums that there is no scaling in the 
coupling constants. There is no good reason why this is OK but this turns out to be 
fine up to n f=s uq. However as was discussed in one of the previous lectures, hidden 
local symmetry with vector manifestation a la Harada and Yamawaki predicts that 
at least the vector meson mass my and the coupling gy run in density proportionally 
to each other as gy/my ~ const, for density n > no- Scalar meson masses could 
in principle scale differently. However the mass of the dilaton scalar Xs connected 
with the trace anomaly of QCD does run like vector mesons in the skyrmion crystal 
matter. This property remains to be confirmed more rigorously. With the BR-scaled 
Landau / functions, the results (second column of (|14p) are considerably improved. 

At n = no 

No BR scaling With BR scaling Exp. 

/3 ~ 18 - 19 MeV 20 - 25 MeV 25 - 26 MeV 
s£ - 0.9 0.7 0.6 -0.7 ^ ' 

K ~ 100 - 136 MeV 140 - 210 MeV 200 - 300 MeV 
8gi - 0.7 0.2 - 0.3 0.23 ±0.07 

The improvement in Sgi is particularly noteworthy. We will return to this quantity 
below in connection with doing Landau fixed point theory in one-go rather than in 
double decimation so far discussed. 

There is a possibility that one can obtain similar improvements by introducing 
three-body and higher-body forces in lieu of BR scaling. To see how 3-body forces 

s N-body forces with N > 3 may also contribute but we will not consider them here. They are 
expected to be suppressed by chiral counting but they could also be suppressed dynamically. We 
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enter in this EFT/RG formulation, it would be necessary to develop a Vi ow k for 
three-body systems. This is because short-range interactions which would be badly 
behaving need to be tamed by an RG procedure as in the case of 2-body force. It's 
not at all obvious how this can be done given the paucity of both experimental and 
theoretical inputs. We know from the C14 dating problem that in certain channels, 
BR scaling and three-body forces do a similar thing and can explain more or less 
similarly the long life-time of the C14. This suggests strongly that they are not 
independent of each oteher. In fact, Song in his thesis showed that BR scaling 
implemented in two-body EFT yields correct thermodynamic consistency conditions 
because what is known as "rearrangement term" involving effective many-body 
forces is implicitly included in the density dependence of the mass and coupling 
constant. 

Nuclear isotopes 



I stable isotopes 

I unstable fluorine isotopes 
I unstable oxygen isotopes 
I neutron halo nuclei 



28 



N 



Fig. 11. Nuclear isotopes to be charted by the RIB machines in the world including the "KoRIA." 
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Fig. 12. Ground state energies of neutron-rich oxygen isotopes. Note the important role of three- 
body forces for neutron number greater than 16. 12SI 



suggest that both BR scaling and 3-body forces should be included together rather than one or 
the other. 
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Let us mention a recent interesting observation that deserves further study. It 
may have a direct relevance to the KoRIA project 0. It has to do with the neutron 
drip line, in particular in the oxygen series 1^21 Shown in Figure [TT] is the chart 
of nuclear isotopes measured and understood as well as those to be measured in 
the future RIB machines. Understanding neutron-rich nuclei will help understand 
the nuclear symmetry energy that is indispensable for the physics of neutron stars 
discussed above. The paper 1=21 addresses the neutron drip line in oxygen isotopes 
and argues that the three-body forces of Fig. 9 are needed to explain the present 
status with the drip line. See Figure [T^J One notes that it is the long-range part 
of the three-body forces that does most of the job, while the shortest-ranged one 
with the coefficient C'e reinforces the tendency. This however seems different from 
the case of C14 dating where the contact term figures more importantly. It will be 
interesting to see what the cause of this difference is. Here is a conjecture: that 
incorporating BR scaling without three-body forces could do the same job in the 
neutron drip line phenomenon. 



3.6. BR scaling in tensor forces 



We know that tensor forces are very important in nuclear physics for various ob- 
servables in nuclear processes. They are considered to be important also beyond the 
nuclear matter density for compact stars, namely in the EOS, e.g., the symmetry 
energy. We know how the tensor forces would scale up to nuclear matter density 
but how the p tensor force scales above the nuclear matter density, in particular, 
from the density at which the half-skyrmion phase (equivalently the "hadronic free- 
dom (HF)") sets in is not known. Note that this is the density regime where the 
"supersoft" symmetry energy turns over. 

Let us see what the "double decimation" strategy says about the scaling be- 
havior. The tensor forces coming from the pion and p exchanges are given in the 
form 



V-= f -^m p ( n .r 2 )(- 
Vj = fy*-m v (T 1 -T 2 )S ia 

47T 



(m p r) 3 (m p r) 2 



3m p r 



exp 



1 



1 



1 



{m^rY (m^r) 2 3m n r 



exp 



(15) 
(16) 



In dense medium, the pion tensor is assumed to remain more or less unchanged but 
the p tensor undergoes the BR scaling. This is of course strictly correct in the chiral 
limit but can only be an approximation when the pion mass is taken into account. 
We will assume that it's OK. In Holt et alP3 where the long-standing C14 dating is 
explained and in Xu and where the "suspersoft" symmetry energy is shown to 
be understandable in terms of BR scaling in the p tensor force, the p mass is taken 
to scale in (|15j) . but the coupling constant f P N is taken non-scaling. 



h This matter was not addressed in detail in this lecture but we think it's worth being added here. 
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Let us see what the current version of BR scaling says. In HLS-type theory, one 

has 

jNp = (1 + K p ) ^ , (17) 

Here K p g P NN is the anomalous tensor pNN coupling and g P NN is related to the 
hidden gauge coupling, so we will simply denote it as g. Thus one is interested in 
the scaling 

f* If „ m *P m N9* ns s 

JNpl JNp ~ — —Z ■ (.1SJ 

Now both empirically and theoretically, the nucleon mass scales up to no as 

m* N /m N « f*/f n , (19) 
and BR scaling says, also up to n , 

m ;/m p » /;//„. (20) 

Thus up to no, 

/^//iV, « ^ « 1 (21) 

where the last equality follows from "double decimation." In 123, this scaling is 
fine since there density only up to no is probed. However in^Sl, this scaling is taken 
beyond no and we will now argue this is not consistent with the modern BR scaling. 

First of all, once one reaches the density corresponding to the half-skyrmion 
density (or HF density), the second equality in (|21[) does not hold. In fact, HLS/VM 
says that it should scale like (qq}*/(qq) which scales as m*/m p = Next, one has 
to have an idea as to how the nucleon mass scales beyond no. In the old version 
of BR scaling, it was taken to scale like $ which means that in the chiral limit it 
should vanish. But this is at variance with the half-skyrmion picture which goes as 
follows. In the skyrmion matter description, the single- nucleon mass in the matter 
must go like 



f> (22) 



where k is an 0(\/ N c ) term consisting in general of a scale-invariant piece (e.g., the 
1/e coefficient of the Skyrme term in the Skyrme Lagranigan) and scale-symmetry- 
breaking (SSB) corrections. If we ignore the SSB corrections, then the nucleon mass 
scales as ()20|) even in the half-skyrmion phase up to the chiral transition point. Now 
we know that in the half-skyrmion phase, /* remains non-zero although the quark 
condensate goes to zero and furthermore it scales little if any. So we may simply let 
the scaling stop at no and take it a constant, 

9 m* N (n ) m p g 
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with 




(my 



(24) 



Taking the empirical value, we will have 



b = m^(no)/mjv~0.7 for n>n 1 / 2 - 



(25) 



This means that once one is in the half-skyrmion phase, the p tensor is suppressed 
by the VM property even though the p mass drops. In fact, a rough estimate shows 
that the p tensor will be almost completely suppressed by n = 3uq, just where the 
"supersoft" symmetry energy drops to negative value. It appears certain that this 
will change the Bao-An Li's scenario on the non-Newtonian gravity. 

3.7. Parity- doubling model for the nucleon 

The observation that the nucleon mass oc /* in the half-skyrmion phase stays non- 
vanishing up to the chiral transition density suggests that the effective chiral La- 
grangian with explicit baryon fields could have a non-z ero chiral invariant mass mo 
as in the parity doubler model of DeTar and KunihirJ^. Work on this picture is 
going on at present, and there is nothing specific to report. 
The general picture, however, is not difficult to figure out. 

For consistency with what we have been doing in this lecture, we need to trans- 
form DeTar-Kuniharu's linear sigma model to non-linear sigma model which is gauge 
equivalent to hidden local symmetry theory that we are working with. To do that, 
define 



in the unitary gauge with £^ R l — £,l — £■ The baryon doubler field VP in the linear 
representation relates to the baryon doubler field Q in the non-linear representation 
through * = with 



where T$ = p^ ® 75. The p's are Pauli matrices acting on the baryon doubler 
(reserving r for the Pauli matrices for isospin) and the 7's a ct on the usual Dirac 
spinor. Thus the SU(2) non-linear version of the model of without explicit chiral 
symmetry breaking is 



£ = 




(26) 



£5 = -r( crl 2 + if ■ 7?r 5 ) = £f 



(27) 



£2 = Qi^Q - 51/ttQQ + g2.UQp3Q 



(28) 



where the vector and axial currents arc defined as 



(29) 
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In (|28[) . the g\ and g 2 terms are the dynamically generated baryon masses and too 
is the chiral invariant mass. One can think of the dynamical mass vanishing as the 
quark condensate vanishes, so in the half-skyrmion phase the baryon mass will be 
given by too. There are attempts in the literature to verify whether this model is 
empirically viable^! So fa r one cannot say. As of today, the reasonable value for 
too is too ~ 0.5 — 0.7 Ge V^^ l Now when applied to dense matter in the sense 
of mean field, if one assumes that the dynamically generated mass goes to zero at 
the density where the half-skyrmion phase sets in, this would give the constant b in 

flUD, 

b ~ 0.5 - 0.7 for n > n 1/2 . (30) 

This is tantalizingly close to what's gotten in the half-skyrmion crystal calculation. 

Although in this lecture we do not go near the chiral transition density for 
which we are not ready yet, there is a very simple and beautiful consequence of the 
scenario given above (i.e., the half-skyrmion or parity-doubled baryon matter) as 
one approaches chiral restoration. It is that the massive baryons will not play any 
role in the thermodynamics of the system. If the nucleon mass remains relatively 
large, say, ~ 700 McV, compared with the meson masses that are dropping a la BR 
scaling, once the density exceeds n^, then the baryons can be integrated out, which 
will thereby endow the parameters of the HLS Lagrangian with density dependence. 
In this case, the properties of the mesons in dense medium approaching the chiral 
transition density n x will satisfy the same formulae as the ones in the temperature 
case except that now T can be replaced by density n. This means that the same 
VM property should hold in density as in the temperature case. This feature will 
play an important role in going toward compact-star regime. 



4. Landau Fermi-Liquid Fixed-Point Theory 

If one accepts t hat nuclear matter is at the Fermi-liquid fixed point as explained, 
e.g., by ShankarfS^l one can get to the answer simply via mean- field theory of EFT. 
To do this, one fixes the Fermi momentum kp, which means fixing the density at 
no. One can then write an effective Lagrangian valid near the Fermi surface as we 
mentioned in connection with Figure 7. Consider having integrated out all degrees 
of freedom other than nucleons and write the EFT action as an expansion in fermion 
field 

S = J d ^f^ (iw - v*k) *P + J UW»H) + 0{{^f) (31) 

with v* = j^j-. Here the expansion is with respect to the Fermi surface at kp 
which is taken to be large, so that A/kp ~ 1/N — > 0. Take the kinetic energy 
term - bilin ear in fermion field - with v* fixed to be scale invariant; then a careful 
analysis shows that the four-Fermi interaction is marginal only when the 4-point 
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vertex U is a constant Uo in momentum space. All momentum-dependent terms @ 
are irrelevant in the sense of RG, so scale to zero by RG flow. The constant Uo is - 
modulo an overall constant - just the Landau / function we used above. The action 
(|31[) is then an EFT for Fermi liquid, with the fixed point describing the nuclear 
matter saturation. All we discussed above will then follow from this Lagrangian. 

Now here is the punch line. The fixed point theory (|31[) is equivalent to a mean- 
field theory a la Walecka of HLS Lagrangian in which fermions are introduced as 
the relevant degrees of freedom §. First gauge-fix the HLS Lagrangian to unitary 
gauge and write Walecka's mean-field Lagrangian as 

C = * [d„ (iff 1 - gw»)\ * + X - (d^4> - m s 2 ) - -F^F^ + ^m^w" + ■ (32) 

We have dropped all other fields (such as pion, p, etc) that can enter in general, 
i.e., in asymmetric nuclear systems. The scalar field <j) in this Lagrangian is not 
the usual scalar denoted as a in the linear sigma model but a scalar \ s that is 
connected to the trace anomaly of QCD (explained by Hyun Kyu Lee and MR 
recently^2l). As it stands, Eq. (|32j) - without filling in the ellipsis - does not give a 
good description of nuclear matter. In fact, it gives too stiff an EOS with K w 500 
MeV. The remedy to this defect is to incorporate BR scaling as we will do below or 
incorporate many-body forces generated by higher dimension operators of the type 

SC = i.9 2 3 + (33) 

where (72,3 are unknown parameters. By adjusting the free parameters of this La- 
grangian, one can obtain, in mean field, an impressive fit of nuclear structure from 
light nuclei to n ucle ar matter. Just to have an idea of how well it works, take a look 
at the reference ^3. We would conjecture that putting BR scaling suitably in (|32|) 
without ([33)) would do equally well. Although this has not been worked out fully, 
the thesis of Song lends a partial support to this conjecture. 

Let us finish this lecture by showing what the mean-field approach with a La- 
grangian of the type (|32| endowed with BR scaling does for the anomalous gyro- 
magnetic ratio considered above. This mat ter is discussed in many places, so we 
won't go into details (which can be found in. I™). Putting in the pion field into (|32|) 
in a chirally invariant way, one gets for the proton 

Sgf = I [^(no)- 1 - 1 - "pionterm" (no)] (34) 

where "pion term" is given as a function of density by chiral symmetry, accurately 
known at low density, and $ is the BR scaling given at n = no- Using $(no) ~ 
■^j" - 1 ~ 0.78 from pionic atoms and "pion term" (no) ~ —0.23, one finds that 

'They contain non-forward scattering amplitude F(<j> ^ 0) in the Landau function. 

J Recall this follows from Weinberg's f theorem as MR discussed in other lectures. The basis of this 

fool 

argument is the work of Matsui J,;> who discussed the equivalence between Landau Fermi-liquid 
theory and Walecka's mean field theory. 
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Sgf « 0.23. This was considered as the first indication that BR scaling is operative 
in nuclei. 

5. Three-Body Repulsion at High Density 

This part was only briefly touched on in the lecture but let us add it. We think it 
will be one of the major issues to resolve if one wants to get a realistic EOS for 
neutron stars. 

The repulsion at short distance which presumably figures importantly at high 
density - and which should be present as is clear from the contact three-body term 
in Figure 9 - must clearly be tamed as was done in the case of 2-body repulsion. At 
present, this has not been done in RG formalisms: There is no 3-body Vi ow k- There 
are indications from HLS/VM in dense medium and also holographic QCD at short 
distances that there should be suppression of such repulsion. In HLS theory in QCD 
as well as in hQCD, the contact three-body force is expected to be given by the 
exchange of vector mesons. At high density, the gauge coupling g will go to zero in 
HLS/VM, so there will be an intrinsic suppression of the repulsion by the vector 
manifestation. In hQCD, as shown by Hashimoto ^, N-body forces are suppressed 
at short distance - and in consequence, at high density - as 0(N c /X N ^ 1 ) where A 
is the 't Hooft constant A = <?y M iV c . Since in hadrons, N c /X ~ 1/10, many-body 
forces will be strongly suppressed. This hQCD argument is indicative of an exact 
many-body solution as it involves a multi-instanton solution that is analytic and 
exact at short distances in the model when the space becomes flat, so in some 
sense indicative of what could come out in realistic theories with a generic 5D YM 
structure (which may be shared by QCD). In going to higher density toward neutron 
stars and ultimately to n x , this is one of the most important issues we will have to 
resolve. 
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